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Description 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates generally to Magnetic Res- 
onance Angiography signal processing and more spe- 
cifically to a method of reducing three-dimensional 
Magnetic Resonance Angiography or stationary tis- 
sue data to a two-dimensional image using statistical 
manipulation, while reducing noise and/or increasing 
contrast. 

2. Description of Related Art 

Known techniques for magnetic resonance an- 
giography include methods that, either by volume or 
sequential-slice acquisition, produce data sets with 
spatial resolution in three dimensions. These techni- 
ques may utilize either time-of-flight effects or mo- 
tion-induced phase shifts to produce flow contrast. 
Since there is at this time no practical real-time three- 
dimensional display technology (e.g., holographic 
video) available, this data is reduced to two spatial di- 
mensions for convenient viewing by the radblogist. 
This has generally been accomplished by generating 
a series of projection Images at sequential view an- 
gles. These Images are then repeatedly displayed in 
rapid sequence on a display screen to give the viewer 
a three-dimensional impression of the vascular anat- 
omy. The projection images are created by ray tracing 
through the volume of data, using some predeter- 
mined method for calculating the projected pixel value 
from each ray. In the case of phase contrast data, it 
is desirable to preserve the quantitative nature of the 
data insofar as is possible. 

The most obvious projection method is to take the 
mean of all volume pixel ("voxel") values along each 
cast ray. This unfortunately yields images with low 
vessel contrast and poor detection of fine features 
and vessels, especially for time-of-flight data that al- 
ways contain signals from unsuppressed stationary 
tissue. This Is because vessels are sparse within the 
data volume and the total noise and/or background 
signal along a ray can easily be comparable to or 
much larger than the total flow signal along that ray. 
To provide good vessel contrast, the method of choice 
has become projection of the largest voxel value 
along each ray. (See Dumoulin, C.L.. Souza, S.P.. 
Walker, M.F.. «nd Wogia, W., "Threo-dlmenslona! 
phase contrast angiography." Magn. Reson. Med.. 
(1989). and Armatur. S.C., Masaryk. T. J.. Modic. 
M.T.: et al., 'SOFT time-of-flight magnetic resonance 
angiography," Dynamic Cardiovascular Imaging, 2. 
170, 1989). This method is generally known as Max- 
imum Intensity Projection (MIP). The MIP method 
uses the maximum voxel value along a projection ray 



and discards the remaining voxel values along that 
ray. Since along each ray which actually intercepts 
one or more vessels the only voxel which contributes 
to the projected pixel is within a vessel, the contrast 

5 is very good. Unfortunately MIP has several draw- 
backs. Since the MIP method retains only one voxel 
along each ray. it does not make efficient use oHhe 
acquired data. This is especially true where a vessel 
extends over many voxels and averaging over the 

10 vessel might be employed to improve the vessel sig- 
nal-to-noise ratio. The MIP method tends to lose con- 
trast for vessels with strong signal near the walls but 
weak signal near the axis of the vessel. This situation 
occurs in phase contrast MR angiography when the 

15 velocity at the center of the vessel exceeds the max- 
imum for which the applied flow encoding produces 
monotonically increasing signal. The MIP method 
also may obscure real pathological features. For ex- 
ample, this will occur when a stenotic lesion distorts 

20 a vessel wall and the lesion appears along the line-of 
sight rather than displayed in profile. There is a need 
for a projection method that retains the maximum 
amount of information possible when transforming 
the acquired data to projection data while introducing 

25 no artifactual features and discarding few actual fea- 
tures. 

OBJECTS OF INVENTION 

30 It is an object of the present invention to construct 

a high-contrast low-noise two-dimensional image 
from three-dimensional magnetic resonance (MR) 
fluid flow data from a volume of a subject. 

It is another object of the present invention to alter 

35 the contrast between structures of a subject imaged 
to make differing structures more visible. 

It is another object of the present invention to 
adapt the MR data to reduce differing types of noise 
inherent in the MR data. 

40 It is another object of the present invention to con- 

struct projection images from three-dimensional MR 
stationary tissue data, thereby producing synthetic 
images similar to conventional x-ray images but with 
alterable contrast between structures. 

45 

SUMMARY OF THE INVENTION 

In this invention, three-dimensional Magnetic 
Resonance fluid flow or stationary tissue measure- 

50 ments of a volume of a subject are processed into 
high-contrast low-noise two-dimensional Images. 
The first step is to acquire by known means a three- 
dimensional data set having a number of entries, 
each entry corres ponding to a Magnetic Resonance 

55 fluid flow or stationary tissue volume element, or vol- 
ume pixel (voxel) value of said subject. The value of 
each pixel in the final two-dimensional image is a cal- 
culated statistic of the array of voxel values intercept- 
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ed by an imaginary ray cast through the three- 
dimensional data along the line of sight to that partic- 
ular pixel at some view angle. An order of the statistic 
is chosen which depends on the general shape of a 
histogram of the data that is to be imaged. A projection 
ray is passed at a predetermined angle through the 
data volume and the three-dimensional data set en- 
tries corresponding to alt intersected voxels are chos- 
en. The average and standard deviations of this set 
of entries are taken. 

The projection average is subtracted from each 
of the entries and the result is then raised to a power 
being equal to that of a statistical order, to obtain the 
weighted voxel differences. 

The sum of the weighted voxel differences is div- 
ided by the number of entries, then divided by a nor- 
malization constant to result in a projection value. 
This procedure Is repeated for all parallel rays cover- 
ing the three-dimensional data set and intersecting all 
pixels in the projection image, and the data set is then 
displayed as a two-dimensional fluid flow image. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The features of the invention believed to be novel 
are setforth with particularity in the appended claims. 
The invention itself, however, both as to organization 
and method of operation, together with further ob- 
jects and advantages thereof, may best be under- 
stood by reference to the following description taken 
In conjunction with the accompanying drawings(s) in 
which: 

Fig. 1 is a block diagram showing the apparatus 
of the present invention that is used for obtaining the 
three-dimensional MR data, and constructing a two- 
dimensional projection image. 

Fig. 2 is a graphical representation of the geom- 
etry employed by the present invention, showing a 
three-dimensional data set volume, a projection ray. 
and a weighted projection value pertaining to the pro- 
jection ray. 

Fig. 3 shows a volume of the subject to be i maged 
that corresponds to the data set volume shown In Fig. 
2. 

Fig. 4 Is a histogram of the three-dimensional MR 
data obtained by a typical phase contrast fluid flow 
image shown in normal amplification and 50 times 
amplifk:ation. 

Fig. 5 is a histogram of a typical "time-of-f light" 
fluid flow image. 

Fig. 6 is a flow chart illustrating use of the present 
invention. 

DETAILED DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

U.S. Patents No. 4,431,968 Edetstein et al. is- 
sued Feb. 14, 1984, and 4.796,635 Dumoulin issued 



Jan. 10, 1989 and assigned to the instant assignee 
disclose an app^ratus such as that shown in Fig. 1 
that is used to create three-dimensional magnetic 
resonance images of fluid flow or stationary tissue in 
5 a subject. Magnet 428 immerses the subject to be Im- 
aged in a static magnetic field. Computer 401 sends 
gradient power pulses via interface 405 to gradient 
power supplies 406, 408, and 410. Gradient power 
supply 406 energizes gradient coil 416. The time- 
to changing current passing through gradient coil 416 in- 
duces a time-changing magnetic field having a gra- 
dient in the x direction. Similarly, current passing 
through coils 418 and 420 produce magnetic fields 
having gradients in the y and z directions respective- 
15 ly. 

Computer 401 sends RF excitation pulse Infor- 
mation to transmitter 402 via interface 405. Transmit- 
ter 402 creates an RF pulse of a predetermined 
strength , duration, and frequency that is amplified via 

20 RF power amplifier 412. The RF pulse is transmitted 
by transmitter coil 424 to the subject. The combina- 
tion of the RF pulse, the static magnetic field, and the 
magnetic field gradients causes magnetic resonance 
of nuclei of the subject that can be localized. 

25 The signals received from the subject are re- 
ceived by a receiver coil 426. then processed by a 
preamplifier 422. receiver 414. and averager 404. 
The signal Is then sent through an interface 405 to 
disk storage 403. The entire set of signals is received 

30 and stored. The signals are reconstructed by comput- 
er 401 into a three-dimensional data set volume 41 as 
shown in Fig. 2. This three-dimensional data set vol- 
ume 41 is composed of voxels, each corresponding to 
fluid flow approximately as it occurs in the original vol- 

35 ume 54 of Fig. 3. 

After reconstructing the three-dimensionatdata 
set 41 . computer 401 represents portions of the data 
on two-dimensional display device 430, via interface 
405. The display device may be a printer, plotter, col- 

40 or- coded color monitor, or gray scale monitor. 

For a more detailed description of operation of 
the apparatus shown in Fig. 1, see U.S. Patent No. 
4,431.968. 

To display an Image on a two-dimensional 
45 screen, a set of projection amplitudes 45 of Fig. 2 
must be taken of the three-dimensional data set vol- 
ume 41 . Casting projection rays 48 through the entire 
three-dimensional data set volume 41 produces a set 
of projection amplitudes 45. Each projection ampli- 
50 tude 45 is created by processing a plurality 43 of data 
set entries 46, each representing a corresponding 
voxel volume of the subject, along projection ray 48. 

The prior art would simply average all data set en- 
tries 46 along projection ray 48 to obtain the projec- 
55 tion amplitude 45 pertaining to that projection ray 46. 

Another method, the Maximum Intensity Projec- 
tion (MIP) method, is applied by simply taking the val- 
ue of the voxel along projection ray 48 having the 
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highest amplitude and using this value as the projec- 
tion amplitude 45 of projection ray 48. 

The present invention employs a method of stat- 
istical calculation of data set entries 46 along projec- 
tion ray 48 permitting an adjustable tradeoff between 
increased contrast (i.e., between small vessels 51 
and background signal as shown in Fig. 3), and an in- 
creased signal to noise ratio. This statistical method 
employs an adjustment parameter k, the order of the 
statistic, which is adjusted to image differing types of 
data. 

Consider a histogram of the voxel values in a 
three-dimensional MR angiogram, as shown in Fig. 4 
for phase contrast data and Fig. 5 for time-of-f light 
data, in each case for one set of rays from a three- 
dimensional data set. In the case of phase contrast 
data there are at least two components to the distrib- 
ution: a large, symmetrical, nearly Gaussian peak 11 
at low voxel values representing noise, and a tail 12 
extending toward high voxel values representing true 
vessel information (possibly with a contribution from 
faint ghosts of the vessels). The noise in peak 11 dom- 
inates the mean of this distribution (its lowest-order 
statistic), with only a small contribution from voxel 
values 12. 

The MIP method takes the maximum voxel value 
(Vmax) along projection ray 48 of Fig. 2. This single 
voxel value Vmax is chosen, and all other voxel values 
are discarded. 

The present invention uses an approach that is 
different from either of the above methods. The pres- 
ent invention uses intermediate-order statistics repre- 
sentative of the distribution of voxel values along 
each projection (e.g.,48 of Fig. 2). Standard deviation 
(or its related measure variance), skewness. kurtosis, 
etc. are defined by: 

where N is the number of samples, V is the 
mean of the distribution, k is the order of the statistic, 
o is the standard deviation of the distribution and V| 
is an individual voxel value Intersected by projection 
ray 48 of Fig. 2. The resulting values are normalized 
to provide the same units as V). In addition, k=2 cor- 
responds to standard deviation, that Is a measure of 
the width of a distribution, k=3 corresponds to skew- 
ness (a measure of the asymmetry of a distribution), 
and k==4 corresponds to kurtosis (a measure of the 
flatness or peakedness of a distribution). Higher or- 
ders of k may be defined, as well as non-integral or- 
ders (e.g.. k=2.7). The MIP method Is roughly equiv- 
alent to k^infinity. 

Each of these statistical measures will indeed de- 
tect vessels and other structures, but the best meas- 



ure is highly data-dependent. In the distribution 
shown in Fig. 4, skewness is likely to be the best 
measure since the high voxel values 12 contribute to 
the asymmetry of the distribution but the low value (or 

s noise) voxel values 11 do not. Using skewness rather 
than the MIP method for the data set yields projection 
images with Improved contrast and signal- to- noise 
ratios for larger vessels without loss of fine detail., , 
In general, lower-order statistics, i.e., lower k val- 

10 ues, yield higher .signal-to-noise ratios, better rendi- 
tion of extended structures such as medium-to-large 
vessels, and better sensitivity to features along the 
projection ray 48 of Fig. 2. They also better preserve 
the quantitative nature of phase contrast data. Hlgh- 

15 er-order statistics yield better contrast against noise 
or stationary-tissue background and better sensitivity 
to small (e.g.. single-voxel) features occunring at the 
right hand side, the tail 12 of Fig. 4. As a general rule, 
one should use the lowest-order statistic that pro- 

20 vides sufficient contrast and background suppres- 
sion. For phase contrast data the best statistic is like- 
ly to be standard deviation or skewness. For tlme-of- 
flight data a higher-order statistic is generally better. 
An example of how the present Invention is to be 

25 applied is shown in Fig. 6. Data is acquired in step 18 
by employing a three-dimensional flow-sensitive MR 
imaging technique. The apparatus of Fig. 1 is enrv 
ployed to collect this data and perform the computa- 
tions. 

30 The characteristics of the distribution histogram 

of the data, similar to Fig. 4 or 5. are constructed in 
step 19. The order of the statistic is determined, which 
reflects the shape of the data in the histogram, (e.g.. 
if the data to be imaged adds to the skewness of the 

35 histogram, the order of the statistic should be "3".) 

The set of data set entries 43. which intersect a 
ray 48 of Fig. 2, is determined. An average set 43 is 
taken in step 21 of Fig. 5. In step 22 the standard de- 
viation of the voxels intersecting ray 48 of Fig. 2 is de- 

40 termlned, except in the case k=2. The standard devia- 
tion obtained is raised to a power equal to the order 
of the statistic minus one to obtain the normalization 
constant at step 23. Each voxel value V| is subtracted 
from the voxel average to create a distribution differ- 

45 ence in step 24. 

The distribution difference Is raised to a power 
equal to the order of the statistic, to obtain a weighted 
voxel difference in step 25. 

It is then determined in step 26 if the voxel value 

50 that was processed was the last voxel intersecting ray 
48 of Fig. 2. If it was not, processing continues along 
path 27, and the next voxel along the projection ray 
48 of Fig. 2 is chosen in step 28. Steps 24 and 25 are 
then repeated for voxels intersecting this projection 

55 ray. 

If the voxel processed at step 25 was the last vox- 
el intersecting ray 48 of Fig. 2, processing passes 
along path 30 to step 32. The weighted voxel differ- 
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ences for voxels intersecting ray 48 of Fig. 2 are sum- 
med in step 32 to obtain a projection sum. The projec- 
tion sum is divided by the number of voxels intersect- 
ing ray 48 of Fig. 2, then (except in the case of k=2)» 
divided by the normalization constant in step 33 to ob- 
tain a weighted projection value for projection rays 45 
and 48 of Fig. 2. (For k=2 proper nonmalization at step 
33 is obtained by taking the square root of the projec- 
tion sum divided by the number of voxels intersecting 
ray48of Fig. 2). 

In step 34 it is determined if the projection ray 
processed is the last projection ray. If it is not the last 
projection ray. the next projection ray is chosen in 
step 36. and processing returns to step 21 to repeat 
the processing for the new projection ray. Processing 
continues as described above until all parallel projec- 
tions through threedlmensional data set volume 41 
(Fig. 2) are processed. The two-dimensional array of 
projection values is then a numerical representation 
of the projection image. 

When all the projections have been processed, 
the processing continues along path 38. The projec- 
tions are then transferred to an electronic two-dimen- 
sional display device (430 of Fig. 1) and are displayed 
In step 39 by computer 401 through the interface 405 
to display device 430 of Fig. 1 device for viewing. 

The methods described herein will be useful in 
obtaining high quality projection angiograms from 
three-dimensional MR flow data. These methods re- 
sult In high-contrast low-noise two-dimensional im- 
ages from collecting three-dimensional magnetic res- 
onance (MR) fluid flow data from a volume of a sub- 
ject 

These methods are also useful when applied to 
MR stationary tissue data to produce projection im- 
ages similar to conventional x-ray images, but with 
contrast between structures that can be altered by 
varying the order of the statistic k. Further, these 
methods are not limited to data obtained by magnetic 
resonance methods, but may be applied to any three- 
dimensional Image data. 

While only certain prefenred features of the inven- 
tion have been illustrated and described herein, many 
modifications and changes will occur to those skilled 
in the art It is, therefore, to be understood that the ap- 
pended claims are intended to cover all such modifi- 
cations and changes. 



Claims 

1. A method of displaying a two-dimensional (2D) 
projection image of desired structures of a sub- 
ject from collected three-dimensional (3D) image 
data of a volume of a subject cnmprlsinn the steps 
of: 

a) acquiring a 3D data set having a plurality of 
entries, each entry being a value correspond- 



ing to at least one physical parameter repre- 
sentative of.a voxel within said volume; 

b) acquiring a histogram of said 3D data set 
and selecting a statistical order k according to 

5 the general shape of the histogram of the 3D 

data; 

c) choosing a viewing angle from which to 
view said volume; 

d) determining a set of intersected voxels 
10 which are intersected by a selected imaginary 

projection ray passing through said volume at 
the viewing angle and detemnining N inter- 
sected data set entries V). conresponding to 
the intersected voxels; 
15 e) calculating an average V and a standard 

deviation a from the N Intersected data set en- 
tries Vi; 

f) calculating a projection value P^ from the N 
intersected data set entries V|. for k not equal 

20 to 2 according to the equation : 

or, for k = 2, according to the equation: 
N 

g) repeating steps "d"-T for a plurality of par- 
35 allel imaginary projection rays at the viewing 

angle covering said volume; ^ 

h) transferring the projection values P^ to a 
two-dimensional display device; 

i) displaying a two-dimensional image created 
40 from the projection values P^ on the two di- 
mensional display device for analysis; 

j) Increase the value of the statistical order k 
and repeat steps "d"-**!** If a greater contrast of 
said structures in the projection image Is de- 
45 sired; and 

k) decrease the value of the statistical order k 
and repeat steps "d"-T if a greater signal-to- 
noise ratio of said structures in the projection 
Image is desired. 

50 

2. The method of displaying a 2D projection image 
as recited in Claim 1 wherein the 3D data set is 
acquired by phase contrast magnetic resonance 
imaging techniques. 

-55 

3. The method of displaying a 2D projection image 
as recited in Claim 1 wherein the 3D data set is 
acquired by ti me- of-f light magnetic resonance 



5 



9 



EP 0 506 302 B1 



10 



flow imaging techniques. 

4. The method of displaying a 2D projection image 
as recited in Claim 1 further including the step of 
varying the order of the statistic for a predeter- 
mined three-dimensional data set to produce im- 
ages with differing contrast in elements in said 
subject imaged. 

5. The method of displaying a 2D projection image 
as recited in Claim 1 wherein the 3D data set is 
acquired by magnetic resonance imaging techni- 
ques used to produce images of substantially sta- 
tionary tissue. 

6. The method of displaying a 2D projection image 
as recited In Claim 1 wherein the 3D data set is 
acquired by tomographic imaging techniques. 

7. The method of displaying a 2D projection image 
as recited in Claim 1 wherein the 3D data set is 
acquired by magnetic resonance imaging techni- 
ques. 



Patentanspruche 

1. Verfahren zur Darstellung eines zweidimensl- 
onalen (2D) Projektlonsbildes von gewunschten 
Strukturen eines Subjektes aus gesammelten 
dreidimensionalen (3D) Bilddaten von einem Vo- 
lumen eines Subjektes, enthaltend die Schritte: 

a) Gewlnnen eines 3D-Datensatzes mit meh- 
reren EintrSgen, wobei jeder Eintrag einen 
Wert hat entsprechend wenigstens einem 
physlkaiischen Parameter, der ein Voxel in 
dem Volumen darstellt, 

b) Gewlnnen eines Histogrammes des 3D-Da- 
tensatzes und Wdhlen einer statistischen 
Ordnung k gemdli der allgemelnen Form des 
Histogrammes der 3D-Daten, 

c) WShten eines Betrachtungswlnket. von 
dem das Volumen betrachtet wird. 

d) Ermittein eines Satzes von geschnittenen 
Voxeln. die durch einen gewahlten imagind- 
ren Projektionsstraht geschnitten sind, der 
bei dem Betrachtungswinkel durch das Volu- 
men hindurchtritt. und Ermittein N geschnitte- 
ner Datensatzeintrdge V| entsprechend den 
geschnittenen Voxeln, 

e) Berechnen eines Mittelwertes V und einer 
Standardabweichung o von den N geschnitte- 
nen Datensatzeintragen V|, 

f) Berechnen eines Projektionswertes P^ aus 
den N geschnittenen Datensatzeintragen V| 
fur k ungleich 2 gem3(l der Glelchung: 



1 N 
Pk = — S (V-V^)^ 

5 

Oder fur k=2 gemaB der Glelchung: 

g) Wiederholen der Schritte "d"-T f Qr mehre- 
re parallel imagindre Projektionsstrahlen bei 

15 dem das Volumen uberdeckenden Betrach- 

tungswinkel, 

h) Obertragen der Projektionswerte Pk zu ei- 
ner zweidimensionalen Darstellungsvorrich- 
tung, 

20 I) Darstelten eines zweidimensionalen Bildes, 

das aus den Projektionswerte n P^ auf der 
zweidimensionalen Darstetlungsvorrichtung 
fur eine Analyse hervorgerufen ist, 
j) Vergrof^ern des Wertes der statistischen 

25 Ordnung k und Wiederholen der Schritte "d"- 

"i". wenn ein groBerer Kontrast der Strukturen 
in dem Projektionsblld gewiinscht wird. und 
k) Verkleinern des Wertes der statistischen 
Ordnung k und Wiederholen der Schritte d"- 

30 "i". wenn ein grSlieres Signa!/Rausch-Ver- 

hSltnls der Strukturen in dem Projektionsblld 
gewunscht wird. 

2. Verfahren zum Darstellen eines 2D-Projektions- 
35 bildes nach Anspruch 1 , wobei der 3D-Datensatz 

durch Phasenkontrast-Magnetresonanz-Bildge- 
bungstechniken gewonnen wird. 

3. Verfahren zum Darstellen eines 2D-Projektions- 
40 bildes nach Anspruch 1 , wobei der 3D-Datensatz 

durch Laufzeit-Magnetresonanzstrfimungs-Blld- 
gebungstechniken gewonnen wird. 

4. Verfahren zum Darstellen eines 2D-Projektions- 
45 bitdes nach Anspruch 1, ferner den Schritt ent- 
haltend, dall die Ordnung der Statistik fur einen 
vorbestlmmten dreidimensionalen Datensatz 
verandert wird, um Bilder mit unterschiedlichem 
Kontrast in Elementen des abgeblldeten Subjek- 

50 tes zu erzeugen. 

5. Verfahren zum Darstellen eines 2D-Projektions- 
bildes nach Anspruch 1 , wobei der 3D-Datensatz 
durch Magnetresonanz-Blldgebungstechniken 

55 gewonnen wird, die angewendet warden, um Bil- 

der mit im wesentlichem stationSrem Gewebe zu 
erzeugen. 
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6. Verfahren zum Darstellen eines 2D-Projektions- 
bildes nach Anspruch 1 . wobei der 3D-Datensat2 
durch tomografische Bildgebungstechniken ge- 
wonnen wird. 

7. Verfahren zum Darstellen eines 2D-Projektlons- 
blides nach Anspruch 1 , wobei der 3D-Datensatz 
durch Magnelresonanz Bildgebungstechniken 
gewonnen wird. 



I 



(V - v-)*^ 



at pour k = 2. d'apr^s I'^quation : 
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Revendications 

1. Proc6d6 pour visualiser une image en projection 

bidimensionnelle (2D) de structures souhalt^es is 
d'un objet ^ partir de donn^es d'image tridimen- 
sionnelles (3D) recueillies concernant un volume 
d'un objet, comprenant les stapes consistant ^ : 

(a) acqu^rir un jeu de donn^es 3D ayant une 
pluratit6 d'entr6es. chaque entr6e 6tant une 20 
valeur qui correspond ^ au moins un param^- 

tre physique reprdsentatif d'un voxel dans le- 
dit volume, 

(b) acqu^rir un histogramme dudit jeu de don- 
n^es 3D et choisir un ordre statistique k en 25 
fonction de la forme g6n6rale de Thistogram- 

me des donn^es 3D, 

(c) choisir un angle d'observation depuls le- 
quel on souhaite observer ledit volume, 

(d) determiner un groupe de voxels Intersec- 30 
t^s qui sont coupes par un rayon imaginaire 

de projection choisi qui traverse ledit volume 
sulvant ledit angle d'observation et determi- 
ner N entr6es V| d'un groupe de donn6es in- 
tersectees qui correspondent aux voxels in- 35 
tersectes, 

(e) calculer la valeur moyenne V et r6cart type 
o des N entries V| du groupe de donn^es in- 
tersectees. 

(f) calculer une valeur de projection P^ k partir 40 
des N entries V| du groupe de donn6es Inter- 
sectees, 

pour k different de 2, d'apr^s ('equation : 
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lite de rayons imagin^ires de projection, pa- 
ralieies, faisant ledit angle d'observation etre- 
couvrant ledit volume, 

(h) transferer les valeurs de projection P^ d un 
dispositif de visualisation en deux dimen- 
sions, 

(i) af f icher sur le dispositif de visualisation en 
deux dimensions une image bidimensionnelle 
creee d partir des valeurs de projection Pk . 
pour analyse, 

(j) augmenter la valeur de I'ordre statistique k 
et repeter les etapes (d) d (i) si on souhaite un 
plus fort contraste desdites structures dans 
I'image de projection, et 
(k) diminuer la valeur de I'ordre statistique k et 
repeter les etapes (d) d (i) si on souhaite un 
plus fort rapport signal/bruit desdites structu- 
res dans rimage de projection. 

2. Precede pour visualiser une image en projection 
2D selon la revendication 1, dans lequel le jeu de 
donnees 3D est acquis par des techniques d'ima- 
gerie par resonance magnetique ^ contraste de 
phase. 

3. Precede pour visualiser une image en projection 
2D seton la revendication 1 . dans lequel le jeu de 
donnees 3D est acquis par des techniques d'ima- 
gerie par resonance magnetique ^ temps de vol. 

4. Precede pour visualiser une image en projection 
2D selon la revendication 1. comprenant en outre 
retape qui consiste k faire verier Tordre statisti- 
que pour un groupe predetermine de donnees tri- 
dimensionnelles afin de produire des images 
avec differents contrastes desdites struclures 
dans ledit objet dont on forme t'image. 

5. Precede pour visualiser une image en projection 
2D selon la revendication 1, dans lequel le jeu de 
donnees 3D est acquis par des techniques d'ima- 
gerie de resonance magnetique utilisees pour 
produire des images de tissus senslblement inrv 
mebiles. 

6. Precede pour visualiser une image en projection 
2D selon la revendication 1 , dans iequet le jeu de 
donnees 3D est acquis par des techniques d'ima- 
gerie par tomographie. 

7. Precede pour visualiser une image en projection 
2D seton la revendication 1, dans lequel lejeu de 
donnees 3D est acquis par des techniques d'ima- 
gerie par resonance magnetique. 
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(g) repeter les etapes (d) e (0 pour une plura- 
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Determine anc average 
vMall voxels intersecting 
oro lection rav 
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calculate standard 
deviation for voxels 
incersectino oroiection 
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^alse standard deviation 
to a power (Deing the 
order of the statistic- i ) 
to get normalization 
constant 
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Subtract next voxel value 
from voxel average to get 
distribution difference 



Raise distribution difference 
to a power (being order of 
the statistic), to get 
weighted voxel difference 
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Sum all weignted voxel 
differences in projection 
get projection sum 
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divide projection sum by 
the number of voxels in :ne 
projection then divide by 
the normalization constant 
to get a weighted 
projection value for one 
projection ray 
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Display all 
projections for 
a 2D image 



Fig. 6 



13 



